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Abstract 

The classical modular equations involve bivariate polynomials that can be seen to be uni- 
variate with coefficients in the modular invariant j. Kiepert found modular equations relating 
some //-quotients and the Weber functions 72 and 73. In the present work, we extend this idea to 
double 7y-quotients and characterize all the parameters leading to this kind of equation. We give 
some properties of these equations, explain how to compute them and give numerical examples. 

1 Introduction 

Let r] denote Dedekind's function. When > 1 is an integer, r/-quotients of the form / = 
Y\d\N vi^ / d'Y'' functions for r'^(A^) when the integer r^'s satisfy some properties known as 
Newman's Lemma |12) . In other words, there exists a bivariate polynomial ^[f]{X,J) such that 
= for all z, where j is the classical modular invariant. 

In some cases, there exist equations of the form <!>[/] (X, 0^,02) where <I>[/](/(z), 73(z), 72(2)) = 
for the Weber function 73, 72. Kiepert was the first to compute modular equations of this type for 
/ = tt)p = r][z / p) / r]{z) for p < 29 (see [9]). Weber cites some examples in [14, §72] and Antoniadis 
[1] extended this to p < 61. 

In the present work, we study such equations for the double r/-quotients rtip^ p^' ^ introduced 
in [5] . We give all parameters {pi ,P2,e) leading to equations in 72 and 73 . 

Section 2 recalls known facts on Weber and tj functions. Section 3 deals with the case of tVp 
where we introduce a faster variant of the classical algorithm to compute the modular equation via 
series expansions. Section 4 proves the necessary results for tDp^^pj, gives algorithms to compute 
the equations in the spirit of Section 3, and we add numerical examples. 

Notations: If u is some function, we will note $[ti](A, J) the corresponding modular equation. If 
u = j{nz), we will note to simplify. 
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2 Preliminaries 



2.1 Properties of the functions 72 and 73 



We will use the traditional notations 



T = 



( 



1 1 
1 



) 



s 



-1 

1 



) 



and use the notation / o M to denote the function z 1— )■ f{Mz). The modular invariant is j{z) = 
j{q) = 1/q + 744 + • • • with q = exp(2i7rz). The classical Weber functions are 



72(g) = iiqf/^ = q~^/^{l + 248g + AUAq^ + 34752^^ + 213126g^ + 0{q^)) G q^^/^{l + 
73(g) = {M - 1728)1/2 = q-^/^{l -492q- 22590g3 + Q {q'>)) G q'^^^l + Z[[q]]). 



If n is an integer, we note (^n = exp(27ri/n). Remember that j is invariant through T and S and 
that 



Moreover, we have [14:\ §55] 

Theorem 2.1 (a) Any function invariant by T and S is a rational function of j. 

(h) Any function f satisfying f oT = —f and f oS = —f is equal to 73 times a rational function 



(c) Any function f satisfying f ° T = C,^ f and f o S = f is equal to 72 times a rational 
function of j. 

(d) Any function f satisfying f oT = —Q^^f and f o S = —f is equal to 737^^ times a rational 
function of j. (Note that — Cs^^ = Cq^^-) 

Let us precise this result in a special case. 

Proposition 2.2 Let T{q) he invariant as in Theorem \2.1[ Suppose that T{q) G for 
an irreducible fraction a/h > with b \ 6. When 6=1, T{q) is a polynomial in j{q). When 6 = 2, 
T{q)/^3{q) is a polynomial in j{q). For 6 = 3, T/72 is a polynomial in j{q) where i = —a mod 3. 
For 6 = 6, T/(7|73) is a polynomial in j{q), where i = —{a + 3)/2 mod 3. In all cases, the 
polynomial in j[q) has integer coefficients. 

Proof: in all cases, the integer i is chosen in such a way that the resulting series T' is invariant 
through S and T, therefore a rational function in j. Noting that T' has integer coefficients, by the 
Hasse principle, so does the polynomial. □ 

From the algorithmic point of view, we have to recognize a polynomial with integer coefficients 
applied to j (q) , given the first terms of the series T{q) ■ Note that we need the order of this series 
to be > 0. We proceed step by step. 

function RecognizePolyInJ(T) 

Input: a series T = Cyq^ + • • • + 0{q^) with integer coefficients, v < and Cy ^ 0. 
Output: a polynomial P{X) of degree —v such that T = P{j{q)). 

1. IZ := T; i := valuation (7^); P := 0; 



72or = C3 72, 720 5' = 72, 
73 7" = -73, 73 o = -73- 



(1) 

(2) 



of j- 
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2. while i < do 

{ at this point TZ = riq^ + • • • + 0{q^) with r , / } 

2.1 P ■.= P + riX-'; 

2.2 n:=n-uj{qy'; 

2.3 i := valuation (7^); 

3. return P. 

Note that we can precompute the powers of j{q) whenever needed, so that each call to the 
function requires O(u^) operations. In large cases, computations can be done using results computed 
modulo small primes and reconstructed via the CRT (as done by Atkin, see [TO]). 



2.2 Formulas for the r^- function 

The following is taken from [6] and will be our main tool in the computations of Section [H 
Theorem 2.3 Let M = ^, | ^ T he normalised such that c > 0, and d > if c = 0. Write 



a = ci 



2^('^') with ci odd; by convention, ci = A(c) = 1 if c = 0. Define 



a \ ^ab+c{d{l-a^)-a)+3ci{a~l)+^X{c){a^-l) 



For X G N write 

ukcl + vkKc = Sk = gcd(a, Kc) = gcd(a, K). 

Then 

where the square root is chosen with positive real part. 

We can decompose the formula into several parts: e(M) = Jac(M)(^24^^^ where we distinguish the 
Jacobi symbol part and the exponent of then, we have the squareroot part Q{M) and the 
r/-part N{M). When dealing with a r/-quotient, the aboves formulas are applied by multiplicativity 
on the different pieces r]{z/d) (see below). 



3 Generalized Weber functions 
3.1 Definition and properties 

Let > 3 be an odd integer. For all factorizations N = ad, let e = gcd(a, d) and consider the 
functions 



(a-i)/2 f 2\ ^q r]{{c + dz)/a) 
\ej r]{z) 



for < c < a with gcd(c, e) = 1. These functions were introduced in [14^ §72]. It is easy to see that 
= z(^-i)/2ry(z/iV)/r?(z) = f(^-i)/2tt)7v(z) where the function von was studied in 0]. 
Weber proves that in all cases, P?^ ^ are roots of a modular equation. In some cases, the results 
are better, for instance: 

Theorem 3.1 //gcd(A, 6) = 1, 12 | c, then the Pcda^2~'^^^ ^^^^ '^'^^ roots of a modular equation. 
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3.2 Computations in the prime order case 

For a prime = p > 3, this setting simphfies to 

Theorem 3.2 The numbers Xc,d,a72(-z)^~^73(-z)^^~"'^^^^ o'^e roots of a modular equation whose coef- 
ficients are rational functions of j{z). In particular, the constant term is (—l)^~^^^'^p. 

Antoniadis extended the results of Kiepert to p < 61 and gave more properties of the polynomials 
[T]. He computed the equation by solving a linear system in the unknown coefficients of the equation, 
using the g-expansion of j{q) and the fact that xo,pi must be a root of the equation. 

A standard approach (already known to Enneper [71 §52]) is to compute the powers sums of the 
roots of the equation, recognize them as polynomials in our variables, and then terminate using the 
classical Newton relations. Inspecting our roots, we see that the g-expansion of a;o,p,i has positive 
order, and all Xi2h,i,p have negative order. So the power sums can be computed using the Xi2h,i,p 
only; we can find formulas for the ^-expansion of Yl^=o ^i2h ip needed. 

A better approach is to look at the reciprocal polynomial, whose roots are the l/xo,p,i and 
l/xi2h,i,p and only the first one contributes to the power sums. Write (p — 1)/12 = e'/5 as an 
irreducible fraction with 6 | 5. Noting that 

p/xo,p,i=q^'~P^/'\l + ---)=q'^'/'il + ---), 

we see that all powers are expressible as functions whose expansions satisfy Proposition 12.21 
The algorithm is: 

1. compute Sk = p/x^ p i and recognize it as a polynomial in the usual variables. 

2. use Newton's formulas. 

3. Remove the powers of p. 

Note that the largest power is (p/xo,p,i)^^^ = q~^^ + • • • ) where the exponent is an 

integer as soon as p > 3. Therefore, we need up to (p^ — 1)/12 terms in the j-series. 
We have 

Si = ll/xo,ii,i = q-^/^ - 2gl/6 _ q7/6 + 2^13/6 ^ 0(^19/6)_ 
Dividing by 7273, we find 

1 + 242g + 0{q'^) 

which must be a polynomial in j{q), hence the constant 1. The other coefficients are given in 
Table m We have replaced 72 (resp. 73) by G2 (resp. G3). The corresponding polynomial is (after 
reductions between variables) 

$[-rti?i](F, G2) = F^^ - G3G2F^^ - 242G^F^° - 19965^3^^ 

-585640^2^*^ + 159440490F^ - 285311670611. 
Taking its reciprocal and removing the spurious powers of p yields: 

_ 990i?6 _^ 440G2F^ + 165^3^2 + 22^^^^ + G3G2F - 11, 
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(11/xo - ■ 



,11,1^ 



Y q-^l'i _ 4g-^/3 + 2gV3 + . . . = ^^(j _ 1244) 

3 _ gg-3/2 ^ 9^-1/2 ^ iQgi/2 + . . . = (^3(j2 _ 1002 J + 59895) 

4 - 8q-^/3 + 20g-^/3 ™ 70^2/3 + . . . = Gf2(J^ - 2488 J2 + 1510268 J - 135655520) 

5 g-25/6 _ iOg-19/6 + 35^-13/6 _ 30^-7/6 _ lOSg-^B + 238^^6 + . . . 

= GzGl{J^ - 2246J2 + 1287749J - 145411750) 
q-^ - I2q-'^ + 54g-3 - SSq'"^ - mq'^ + 540 - 418g + • • • 

= J5 - 3732 J4 + 4586706J3 - 2059075976 + 253478654715 J - 2067305393340 

7 g-35/6 _ 14g-29/6 _^ 77^-23/6 _ 182g-17/6 + 924^-^/6 _ 1547^1/6 _^ . . . 

= G^G2{J'' - 34907"^ + 4063139 J3 - 1796527998 + 247854700555 J - 4740750382830) 

8 g-20/3 _ i6g-i7/3 + i04g-i4/3 _ 320^-11/3 + 260^"^/^ + 1248g-^/3 - 3712^^2/3 
+1664^V3 + . . . = Gf2(j6 _ 4976 J5 + 9210680 J'' - 7786404608 + 2955697453292 

-418137392559040J + 12629117378938720) 

g-15/2 _ i8g-i3/2 i35g-ii/2 _ 510^-9/2 ^ 765^-7/2 _^ 1242^-5/2 _ 7038g-3/2 + 8280g-^/2 

+9180g^/2 + • • • = G2.{J'^ - AlUJ^ + 8386065J5 - 6877048710J^ + 2611195915626J3 
-398512009001700 J2 + 16457557949779815 J - 41283301866181650) 

10 - 20g-22/3 + 170g-^9/^ - 760g-^6/3 + 1615g-^^/^ + 476^"^°/^ - 11210g-^/3 

+22440g-''/3 + l615g-^/3 _ 64600^2/3 + . . . 

= G2{J^ - 6220 + 15382190 J'^ - 19242776200 + 12809764457825 J'' 
-4368737795118764J3 + 6696193526329257507^ 
-33921007872189625000J + 233702090524237500000) 

11 g-55/6 _ 22g-49/6 + 209g"'^3/6 _ io78g-37/6 _^ 2926^-3^/^ - 1672^-25/6 _ 15169^-^9/6 
+47234g-i3/6 _ 3i350g-7/6 _ 107426^-^/6 + 218680^^/6 + • • • 
G^GliJ^ - 5978J^ + 14256527 J6 - 173121086707^ + 1 1327366012605 
-3889904574252522J=^ + 631 138185556080950 
-38141443583282670180J + 473098671409604281800) 

12 - 24g-9 + 252g-^ - 1472g-^ + 4830g-6 - 6048^-^ - 16744g-'' + 84480g-3 
-113643^-2 - 115920g-i + 534612 - 370920^ + • • • 

jio _ 7464J9 + 23101236 - 38353325536J^ + 36913772324730J6 - 20784851556729552 
+6580486714450069928 J4 - 1063011399511905159360J3 + 72005127765018136775955J2 
-1322204967509387392211000J + 1424583710586688670191932 



Table 1: Computations for p = 11. 
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already computed by Weber. 

Note that one drawback of the approach is the large degree and sizes of the coefficients before 
reduction via Newton formulas. However, if computations are performed using CRT primes, this 
is not a problem, since we compute the final polynomial modulo the primes. 

The smallest cases are 

$[tT.i] (X, G2) = X^ + WX^ - G2X + 5, 

$[-tt)?](X, G3) = X^ + UX^ + 63X^ + 70^2 + G^X - 7, 
$ [lt)?3] (X, J) = + 26 X^^ + 325 _^ 2548 X^^ + 13832 + 54340 + 157118 X® + 333580 X'^ 
+509366 X^ + 534820 X^ + 354536 X^ + 124852 X^ + 15145 X"^ + (746 - J) X + 13. 

Remark. We concentrated here on the prime index case. The same work can be done on composite 
ones. Note also that we could use resultants for that task, noting the following. Suppose p is prime 
and M is an integer prime to p; write N = pM. Write 

tt'pM(^) = {r0p{z)X0M{z/p)Y ■ 

On the other hand: 

^tv;^]iw;^iz),jiz)) = 0, 

^tv%]{w%iz/p),jiz/p)) = 0, 
%{j{z),j{z/p)) = 0. 

Writing Z = Wpj^{z), X = tVp{z), Y = xvm{z/p)-, the different quantities are related via the 
algebraic equations: 

Z = X'Y\ 

$[ra;^](x^\j) = o, 
$[tt)^](y^^J') = o, 
<i>p( J, /) = 0, 

and the variables can be eliminated via resultants to get a modular equation in Z and J, that needs 
to be factored to get the correct polynomial. 



4 Double ?7-quotients 

4.1 Definition and statement of tiie result 

For primes pi and P2, let 

where s = gcd(24 (pi^-i)(p2-i)) smallest integer such that sr is an integer, where r = {pi — 

l){p2 — l)/24. Note that s \ 24; and s \ 6 when pi and p2 are odd primes. It is shown in [6] that 
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the function tt)p^ is a function on T^{piP2)', properties of the classical modular equation are also 
given. 

We can now state the result that we will prove in this Section. 

Theorem 4.1 Letpi,p2 be two primes, N = piP2, s = 24/ gcd(24, {pi — l){p2 — l)), s a divisor 
of s and 6 = s/e. If N = 1 mod 6 and the parameters are chosen in TablelM, then there exists a 
modular equation ^[[—lY^^tVp^ p^] whose coefficients are rational functions in 73, 72. 



Pi 


P2 


s 


e 


6 


2 


2 


24 


8 


3 


2 


5 mod 12 


6 


2 


3 


2 


11 mod 12 


12 


4 


3 


3 


3 


6 


3 


2 


3 


7 mod 12 


2 


1 


2 


3 


11 mod 12 


6 


3 


2 


5 mod 12 


5 mod 12 


3 


1 


3 


5 mod 12 


11 mod 12 


3 


1 


3 


7 mod 12 


7 mod 12 


2 


1 


2 


7 mod 12 


11 mod 12 


2 


1 


2 


11 mod 12 


11 mod 12 


6 


1 


6 



Table 2: Values of pi and p2 leading to a modular equation ^[(— 1) "•"^rOp^ p^]- 

The following Lemma is used in the Theorem. 

Lemma 4.1 Let 6 G {2, 3, 6} be as above and suppose N = pip2 = 1 mod 5. Then pi = —1 mod 5. 

Proof: For 5 = 2, N = 1 mod 2 gives the answer. When 3 \ 6, we cannot have pi = 3 since 
= 1 mod 6. For 5 to be equal to 3 (resp. 6), surely we cannot have pi = l mod 3 (resp. 6). This 
leaves pi = —1 mod 3 (resp. 6). □ 

The proof of the Theorem will take use several intermediate results that we will present in as 
much a compact way as possible. When pi ^ P2, we will make the convention that pi is odd (so 
that we may have p2 = 2). Moreover, we let u and v be two integers such that upi + vp2 = 1. To 
simplify the proofs, we will be mostly looking at properties using p2, this case being complicated 
when p2 = 2. Reciprocally, using pi and p2 supposes that pi ^ P2- The results and proofs are of 
course symmetrical by exchanging pi and p2- In case of equality, we will write pi = P2 = P- 

4.2 The conjugates of rOpj 

In [6] are given the conjugates of rOp^p^ (with some minor typos). Here, we need precise the 
expansions of tt5pi,p2- Iii view of Theorem 12.31 the value of tDp^^pj o M can be composed as 

»Pi,P2 °M = Jac(M)C24^^Q(M)AA(M) 

where the first part cumulates Jacobi symbols, the second the exponents of (^24, the third one is 
the product of the squareroots and the last one the rj quotient. To ease notations, we also put 
(j) = C24'' = C24^ ^^^^^ ■ We use the notations and philosophy of computations from [6]. 
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Proposition 4.2 Let pi and p2 be two primes. In all cases, we have the following conjugate 

functions: 





M 




tt)pi,p2 O M 


ord I 




(J 







T / — vr 
PlP2 




' 


o) 




—r 1 



The remaining pi + p2 conjugates are: 

M tOp^^p^oM ord I 



M- 



1,0 



WP2 -WPl 
1 1 



2,0 



1 1 



where in the case of C2,u, we set v = —{npi) ^ modp2>^^2 = (1 + Pi^^i^) /P2 for fj, ^ (equivalently 
u ^ 0; fi = corresponds to v = 0). When v > 0, we get 

^^^"^"l (3pi+2)i^ ^fp2 = 2. 



Moreover 



Also, 



3 in\ ) ^^(^"2 + l)+u-l ifp2^2, 

h{0)-^ (3n+2)(.-l) ^^^^^2. 



£2 



/'PI) 

P2'' 

1 ifp2 = 2. 

When pi = P2 = p, we must consider the p — 1 following conjugate functions: 

matrix ^p,p ° ^ ord / 

where 1 = —^v + vp, e{v) = (■^) if p odd (resp. 1 when p = 2) and 

pv (1 — p?') + (— 3p + 2 + ti) ^ — 3 + 3p if p is odd, 



ifp = 2. 



Proof: the cases of the Ay matrices and of B are treated without difficulty, as in [6] . The value of 
Q{M) is 1, unless we are dealing with the case pi = P2= P- The computations for the C matrices 
involve non-zero exponents for ^24- 
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Case pi y^p2- 

In the same lines as in [6], we first prove the result for C2^u when ly > 0. Iterate over 1 < /U < p2 
and define = — modp2 S {Ij • • • )P2 — 1}, V2 = {1 + ^vpi)/p2- Note that i/ i— )■ /i is an 
involution and the corresponding U2's are equal. Moreover, iterating over 1 < jj, < p2 is the same 
as iterating over 1 < v < P2- We find 

^''^ r]{z)r]{pi{z + u)/p2) vi^^ + ^)viPi{z + i')/P2) 

(a) Assume first p2 / 2. We compute Jac(M2,^) = (|^), and the total exponent of C24 is 

1^(1 - pi) + £{M2^^) = {pi-l)[p2iypifJ.'^ +P2 1^ + 2^1 P2- fJ-V2-i' - fj) 

= {pi - l){{p2 - l)(i^/uV + i^ + 2fi) + n{l -V2 + p.ypi)) 

= {pi-l){p2-l){l^fJ.'^Pl + l^ + 2fI + fJ.V2) 
= {pi-l)ip2-l)iKiP2 + l)v2 + l) + l^) 

where we have used P2V2 = 1 + fJ.i'Pi twice. 

When p2 = 2, we find Jac(M2,^) = 1 and the total exponent of C24 is 

u{l - pi) + £{M2,f,) = {pi - 1) (3pi^2(^ ^ 1) ^ ^(3^ _ 1) ^ ^)/2. 
Since v is odd, n = 1 and the exponent reduces to 

(pi-l)(3pi + 2)^. 

(b) For C2,o, 

7/(2: + i)7/(piz/p2) 'n{z)'n{piz/p2)' 

Assume first p2 7^ 2. Then Jac(M2,o) = (||) and the exponent of C24 is 

Pi-1 + £{M2,q) = -ipi-l)i{p2-l){piu^ -2u + l)+u{piu + v-l)) 
= —2Ar{piv? — 2u + 1 — uv) 
= 24r{uv{p2 + 1) + u - 1). 

When p2 = 2, we find Jac(M2,o) = 1 and the total exponent of C24 is 

Case pi = P2 = P- 
In all cases: 

r]{pz)'^ 



rj{z)r]{z + i^/p) 
where 1 = —fii/ + vp. 

When p 7^ 2, we find Jac(M^) = and the exponent given by 0{i'). When p = 2, Jac(M^) = 1 
and the exponent given by v — 1 = 0. □ 
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4.3 Action of T and S 

This section is devoted to the proofs of the actions of T and S on our basic functions as stated in 
the following two propositions. 

Proposition 4.3 (i) BoT = (t)-^B. 

(a) For < 1^ < N - 1, we have A^oT = A^+i; Ajy^i oT = (/)"Mo- 

(in) ForO<u <p2-l, Cs,^ o T = /2(^)-e2(^+i)c<2^^^^ . c^^^^^^ o T = /2(P2-i)-e2(o)+i^2,o. 
(iv) Forl<u <p, C^oT = Clr'^Cu. 



Proof of Proposition | 
(i), (ii) and (iv) are direct applications of Theorem 12.31 

(iii) For < < P2 - 1, one has oT = /2M-e2(z^+i)(^2,i.+i. For u = p2 - I: 

n ^(^)^(^^i(^ + P2)) _ ^, r? + l) ?? + iV) 
V{z + P2)V (^^^^) viz +P2)v{^+Pi 



Proposition 4.4 For all primes pi and p2, one has 

(i) (Ao,i?)o5= (5,^0). 

(ii) When < < p\P2 and gcd{iy,pip2) = 1: 

A^oS = /a^'^^A 

where 1 = —uv + ^12(^1^2) and 



-Lou^ - 2u + uj + 3 + UV12 i/p2 7^2, 
oj + i^{vi2{l-2pi) + 2) ifp2 = 2. 



Suppose from now on that pi ^ P2- The following hold: 
(iii) When < = pip < P1P2: 

A,oS = <i)'^^P^C2,^ 

where 1 = —wv + wp2, and 



94{p) -- 

(iv) When < < p2, use fi = —l/{vpi) mod p2 and 



e2{Tu)+ p{w{p2 + l) + l) + w ifP2^2, 
+ 32i±ip2 + p(3^ _ 2) + ifp2 = 2. 



C2,u o S — A^p-^ . 

(v) We have 



^fP2^2, 
ifp2 = 2. 
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Proof: 

(i) We first get: 



^Pi,P2 ° {T" oS) = Wp,,p^ of^"^ 



and the case u = yields immediately Aqo S = B. On the other hand, we also have the reassuring 
result that 

BoS = K>pi^p^{-N/z) = tt>pi,p2{z) = Ao{z). 
(ii) When gcd(z/,pip2) = 1, we write 1 = —lov + vi2{piP2), and find 



A/ (,M2,;u O = ■ — = (,24 A 

ri{z)ri{{z +uj)/pi/p2) 



When p2 7^ 2, Jac(M2,yit o S) = 1 and the total exponent of (^24 is 

oj + £{M2,^,oS) = -a;i/2(24r-l) + i/(-48r + l + ?;i2(l-pi-p2)) + 24r(a; + 3) 
= 24r(-a;z^2 -2v + uj + Z) + v{ujv + 1 + vi2{l - Pi - P2)) 
= 2Ar{-ujv^ -2V + 0J + Z + VV12) 
= 24r^3(i/). 

When p2 = 2, we also have Jac(M2,y:i o S) = 1 and the exponent of C24 becomes 

{pi~l){uj + u{vi2{l-2pi) + 2)). 

The same type of computations show that the results also holds for p\= P2 = 2. 
(iii) Suppose now that v = ppi, 1 < p < P2- We write 1 = —wiy + wp2- In all cases 



N{M2,t, o S) 



V{z)r]{{pi{z + zu))/p2 



ri{z + w)r]{{pi {z + m))/p2) 
= (2!'^^^ {r''^^^e2C2,^) . 

Assume p2 / 2. Wc get Jac(Af2,^ o S) = (||). The partial exponent is given by 

-w{pi-l)+S{M2,^,oS) = ipi-l){p2Piwp'^ + i-w + 2p2-l)p + ip2-l)Tu) 

= {pi - I){p2p(pip1^) + {-W + 2p2 - l)p + {P2 - 

= (pi - I)ip2p{wp2 - 1) + {~W + 2p2 - + ip2 - i)-!^) 

= (Pl - 1)(P2 - l)ip{w{p2 + 1) + 1)+W) 

yielding the final result. 

When p2 = 2, we find Jac(M2,^ 5) = ! and 

-w{pi - 1) + S{M2,^, o S) = [pi - 1) (p{3w - 2) + S^i^p^ + ^\ _ 
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(iv) For 1 < 1/ < p2, we compute fj, = —\/{ypi) modp2 and 



(v) In all cases, we compute 



Jac(M o S)Af[M oS)= — — — j— = — C24 "TTT T^- 

\pij ri{z - l)r]{p2z/pi) \pij ri{z)ri{p2z/pi) 

When p2 7^ 2, this yields 

Jac(M o S)J\f{M oS) = C24"^'(^"^i^°)Ci,o. 

The exponent of (^24 is 

1-P2 + S{M oS) = {p2 - l){piP2v'^ + {u + I - 2pi)v + pi - 1) 
= {Pi-l){p2-l){-uvipi + l)-v + l) 
= 2ir{-uv{pi + 1) -v + \) 
= -24r6li(0) 

so that C2,o oS = </'-^^i(°)Ci,o. 

When p2 = 2, the exponent of C24 becomes 

-1 + £{M o S) = {p, - 1) + ^ 

so that the final answer is 



Proposition 4.5 We suppose that pi = P2 = P- Then 

(i) When v = pp, 1 < p < p, set 1 = —wp + wp. Then Ai,o S = C^. 

(a) For all p, and all \ < u < p, one has Ci, o S = A^p where p, = mod p. 

Proof: 

(i) When p^2: 



\pj r]{z)r]{z + w/p) 

p2pCT+(-3p+2+«))p+pro-3+3p 



= V^[-y) {yV^^^XM^'^^c^) C24 

^— p^pw+{— 3 p+2+w) p+p-uj— 3+3 p—6{w) ^ 

— ^24 

using Q{w) =pw{l- p^) + (-3p + 2 + w) p-^ + ^p. 
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When p = 2, p = 1 implying tu = w = 1 and 

A2oS = V2C24 / X . . 77^ = C24 = ^-'l- 



(ii) In all cases, we get 



ri{z + pp)r]{{z + pfj,)/p'^) 



4.4 Finding invariant functions 

The idea is simple. Using the explicit actions given above, we need to find suitable modifications of 
the functions B, A^, Ci^^, C2^u such that the action of T and S on any power sum coincides with 
the action on 72, 73 or the product 7273, as given in Section [^31 

Note that oT = (^^^^'"^i?^. Write re = t/6 and remark that this fraction is irreducible (s 
being prime to t implies 6 is). This leads to set x = 4>~'^ = (24'^^'^ ~ C^"*; & primitive 5-th root of 
unity. 

The aim of this Section is to prove the following Theorem from which Theorem 14.11 will follow. 



Theorem 4.6 Assume we are in the conditions of Theorem \4.1\ Define the functions 
— X a — a , — X '-'2,1/ — A "-^2,1/) '^u — X 



where ij, = —l/v modp and 

( 1 ifS = 2, 
ao = < if S = 3, 
[ 3 tfS = 6, 

making x"" = (—1)'^+-'^ = x "° = X 

Then, for all integers k, the quantity 

N-l Pi-1 P2-I 

■Sk = B'^ ^'J^ + X] ^[ J' + ^ C2^J' = B'^ + SA,k + Scx,k + Sc2,k 

jy=0 I/=0 V=Q 

satisfies Sk o (T, S) = (x^ x""'')5fc. 
With these notations, we have 

Proposition 4.7 The following hold: 

(a) B'oT = xB'. 

(b) {A',},oT = {xK},. 

(c) {Cl,}.oT = {xCl,}u. 

(d) For allv, C^oT = xC^. 

Proof: (a) and (c) follow easily from Proposition 14.31 

(b) We first obtain A^j^_^ oT = x^o- explain how the choice A'^^ comes from. For some 

function a to be precised later, let us put A'^ = x'^^'^^-^t-, for which 
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We must find a s.t. 

a{u) -a{v + l) = l mod (5, < z/ < TV - 1, 

and 

a{N - 1) - a(0) + 1 = 1 mod 6. 

The first set of equations gives us a{v) = a{0)—iy mod S and the second a(0) — (A^— 1) = a{0) mod 5, 

which is possible only when = 1 mod 5. Setting = a(0) yields the result. 

(d) Proposition O gives us C^, o T = C2a^~^^CI. A glance at Table [2] shows that p'^ - 1 = 

mod (24/e), which implies 2e{p — 1) = —{p — l)^e mod 24 and therefore C24^^ = X- ^ 
The actual value of a is in fact dictated by the other invariance properties that follow. 

Remark. This proposition shows at the same time that we cannot expect some nice T-action when 
N ^1 mod 6. 

Let us turn our attention to the S'-action on our candidate functions, using the notations of 
Proposition 14.41 



Proposition 4.8 (i) {B', A'q) oS = x'^'C^o, B'). 
(ii) When gcd{v,pip2) = 1, A'^ o 5 = x"°^L- 
(in) For v = pip, A'^o S = x""C£„. 

(iv) For I < V < p2, p = —l/{vp{) mod p2 and €'2^0 S = x"°^^pi • 

(vi) For V = pp, I < p < p, set 1 = —wp + wp. For all p, A'^o S = x^^'C'^. 

(vii) For 1 <u <p, setting p = modp, we have C^o S = x'^'^A'^p- 

Proof: 

(i) We have B' o S = x~"°A'q; A'qO S = x"°B' and the result follows from x~°'° = XQ^o- 

(ii) Proposition 14.41 can be rewritten 

Kos = x''-''-''^''^Al 

and we simplify the exponent using 1 = —cou + V12 mod 5, which leads to: 



A'oS = A' 



When p2 7^ 2, we use X ^ = X"°- The p2 = 2 can occur only for (5 = 3, in which case pi = —1 mod 3 
and the exponent is x is 0. 

(iii) For v = ppi: we use 1 = —wu + wp2 to get 

AloS = x-''^'^Cl^ 

or 

x-'^'+'K ° s = x-'^^'h^-'^^^^c'^^^ 

and we need simplify: 

-Ba{p) - 62(117) + aQ + w - V. 
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Using the definition of 04, we get 

A'os = a. 



2,t^| ^ao-p((3?^l±l)p+pi+3«,-2) .f^^^2. 



and we conclude using pi = —1 mod (5. 

(iv) For 1 < 1/ < p2, we compute fj, = —\/{ypi) mod^?2 and 

Simplifying the exponent gives 



CL„oS = A', 



X 



-a()+At{(p2+l)«J2+Pl+l) jf p2 ^ 2, 



2,z/ ^/ipi 1 ^-Qo+4pi+l iip2 = 2. 

where for p2 = 2, we used v = ^=1. We conclude as in (iii). 
(v) When p2 7^ 2, we start from 

from which 

or 

and the exponent is 

uv{pi +P2 + 2) +U + V -2. 

This quantity is = u + — 2 mod 5 since = — 1 mod 5. Moreover 1 = pi{u + v) mod 5 and finally 
the exponent is —3 mod 5. 
When p2 = 2, we have 



X 

or 



(pj+Pl-3)u^-3 

C^,oo5 = X ^ C-i^o, 

and this is since this can only happen when (5 = 3. 

(vi) Since 1 = —wp + wp, we can write 

and the result comes from the definition of n. 

(vii) C^oS = A^^or 

c:,os = x'^('')x-"°+'^^^;.p, 

and we conclude as in (vi). □ 
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4.5 Properties of the modulEir equation 

Prom the preceding sections, we see that 

N—l P2-1 
i/=0 1^=0 u=0 

is a modular equation whose coefficients can be expressed in terms of j, 72 or 73 depending on the 
value of 6. Before doing this, we may express these coefficients as Puiseux series. 

Proposition 4.9 With the usual notations: 

(a) the coefficient of smallest order of ^ is q~'^^'^; 
(h) the trace has order re; 

(c) when pi / p2, ^{0) = 1/ 

(d) when pi = P2 = P, ^(0) = (— when p is odd and 2^ when p = 2. 

Proof: (a) the coefficient of smallest order comes from the coefficient of pi'W-^-'^ which has the 
order of B' n^=o^ K or 

N-l 

E—re 
— — = —2re. 
N 

When pi p2, — N — 1 = P1+P2; when pi = p2 = p, this is p — 1. Note that all other terms 

have orders stricly less than this bovmd. 

As an example, when s = e, the degree of the equation in J is 2rs and the corresponding term 

is j2rspi^{N)-N-l_ 

Moreover 

N-l /N-1 \ 

^' n< = n x""-^ cr^^^-^)/\-^-(i + • • • ) 

_ Nao-N{N-l)/2^-2ireN{N-l)/2 2re(-, , n 

— ^,ao-NiN-l}/2^-24reN(N-l)/2 2re(-, , 

— A ',24Ar y V-*- "I J' 

using N = 1 mod 6. When N is odd, this reduces to 

^ao-iV(JV-l)/2+(JV-l)/2^-2r-e^j + •••)= ;^°o-(^-l)V2^-2r-e(i + •••). 

(b) The dominant term in the sum of the conjugates is that of B', namely q~^^. 

(c) Por pi^P2-- 

n ^2,'^ = n x^'^''^"''x"^'^''^£2Cpr9''''^^'(i +•••)= x~^^^^~^^^^4^XpT^^^~^^^^Q^'^i^ +•••)• 

Multiplying all together, we find the norm to be of valuation 0, hence a constant 

,Q _ ^ao-Ar(Ar-l)/2A-24reAr(Ar-l)/2 pi(pi_i)/2 «ie>repi(pi-l)/2 -p2(p2-l)/2_P2e>rep2b^ 
^ — X S24Ar A =•! Spi A <^2 Sp2 
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_ /^i^2xe^ao-Ar(Ar-l)/2-pi(pi-l)/2-p2(P2-l)/2A-24reAr((iV-l)/2-(pi+p2-^^ 
— 1^1 ^2 ) X S24iV 

When p2 = 2 (with pi odd), we have 6 = 3 always, meaning ao = and N = 1 mod 3. 
Therefore, noting that e is always even: 



since pi = —1 mod 3. 

When p2 7^ 2, both pi being odd, we may use the quadratic reciprocity law to find 

^ ^ ^_j^e(pi-l)(p2-l)/4^ao-Ar(Ar-l)/2-pi(pi-l)/2-p2(p2-l)/2^^24r-eAr((Ar-l)/2-(pi+p2-2)/2) 

Since pi + p2 — 2 is even, we obtain 

„ _ >3(24re) ao-Ar(Af-l)/2-pi(pi-l)/2-p2(P2-l)/2+(Af-l)/2-(pi+p2-2)/2 
~ S24 A 

^ ao-3-Ar{Ar-l)/2-pi(pi-l)/2-p2(p2-l)/2+(Ar-l)/2-(pi+p2-2)/2 

A. 

_ ^-Ar(Ar-l)/2-pi(pi-l)/2-p2(P2-l)/2+(Ar-l)/2-(pi+p2-2)/2 
A ) 

and by inspection, this is always 1. 
(d) When pi = P2 = P, we get 

n = n x'^^''y/^e(.)^c2i^'^^c2-4r9^^-'^/''(i + • • • ) 

i/=l u=l 



e 



— p-l, 



\i/=i / 

The quantity H^Zi is 1 for p = 2; when p is odd 



TT < ^ ({-\Y-\v-\)\\ f-l 

11 = 



p J \p 



using Wilson's theorem. 

When p = 2, e = 8, (5 = 3, we find 

,9 _ ^ao-N{N-l)/2^-2ireN{N-l)/2 p{p-i)/2 e{p-l)/2^eYZ':X e(u) .-ep{j>-\)/2 
" — X S24Ar A P S24 '524p 

_ .,^0-6/--8-4-3/2 i„4/~8 _ o4 
— A S96 A ^ S48 ~ ^ • 

When p is odd 

,Q _ (-'^\ ''^ao-iV(JV-l)/2+p(p-l)/2 p(p-l)/2^-24reiV(JV-l)/2^e(-(p-l)/2+E^li ©H) 
— I ~ I A P S24Ar S24 

Now: 

p— 1 1 p~i 

= ^pz/(l-^2)_^(_3p + 2 + ?;)/x-3 + 3p= (p-3p + 2)5o + 3(p-l)2 + ^-pz//x2 + 'y/x 

17 



where 5o = X]!^=i ^- Using 1 = —fiu + vp, the sum becomes 

p— 1 p— 1 

— ^p) + = /u(p — i^p^ + v) = pSo mod 24 
j/=i u=i 

in ah cases: when p > 3, = 1 mod 24; when p = 3, u = 1 (resp. u = —1) leads to /i = —1, u = 
(resp. = 1, f = 0). Therefore, the exponent of C24 is 

= e{-{p - l)/2 + {-p + 2)So + 3(p - 1)2) = -e(p - 7)(p - 1)V2 mod 24, 

so that 

'^^''pe(p-l)/2^"0+{A^-l)/2-Af(iV-l)/2+p(p-l)/2+(p-7)/2 ^ / ^ype(p-l)/2^ao-(p4-3p2_8)/2_ 



For instance, when p = 3, e = 3, (5 = 2, we find •& = (— 1)3^(— 1)^^^^ = —3^. More generally, as 
soon as p > 3, 

since = 1 mod 24 and the fact already used that x"" = X~'^- ^ 

4.6 Computing the modular equations using series expansions 

There a variety of methods to compute the modular equations. For large computations, it is possible 
to use suitably modified versions of [3] or [2]. Also, we can use resultants in the same spirit as in 
the remark at the end of Section[3l noting that ^p^^^p^ = i^piiz)/^piiz/p2)y ■ 

Here, we content ourselves to use series expansions and nice formulas that can help us for small 
cases. Also, this will add new properties to our equations. 

Looking carefully at the expression for 5^ , we see that the terms in Ci , C2 or C cannot contribute 
to the modular equation, since they have positive order. Therefore, we need only consider the 
expansions of B''' and SA,k- Doing this, we see that the useful terms for SA,k fii^e for j < —ktN'/N. 
Since B'^ = q~^'^^{l + - ■ ■ ) and 1 < A; < ^(A^). We need at least reil^{N) terms in the last coefficient. 
Since B' is the product and quotient of very sparse series, it might be worthwhile to compute its 
powers by successive applications of special routines handling this kind of computations. It is 
possible to compute nice formulas for the 5^,*;; in the spirit of the ones to come, but we do not 
need them. 

A second algorithm consists in grouping 

<^{F) = Pb{F)Pa{F)PcAF)Pc,{F) 

and to compute Pa (resp. Pq^ and PC2) via its power sums that are given in the preceding 
propositions. 

Inspired by the approach of section 13. 2| the third algorithm uses the reciprocal polynomial, 
whose powers sums will depend on the C{ ^ and C'2 only: 

P2-1 1 pi-i 



/ ^ (^1 k ~^ / J f^i k ' 
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which is a process involving pi + P2- We will prove two useful results (propositions HTTOl and H.lll 
below) to help us compute these quantities. 

Proposition 4.10 For all integers k ^ 0, 

j>ktp'^/p2 

where {p2 + l)/5 = p'2 and the Ck^i are explicitly given in the proof. 
Proof: put w = q^^P'^, C = ("pj and write 

with Ci2{q) = 1 + ■ ■ ■ G 1j[[q]] (which is symmetrical in pi and P2), from which 

= x''^''^-'x-''^''^e-2iwCY'Ci2{wCr. 

and 

P2-1 

c ^ke„,,kTe \ ^ ^—kv/-kreup /'»,,A^^^efc 

Writing Ci2{wY'^ = Yl'iZo'^k,iW^ (remark this is valid irrespective of the sign of k), the inner sum 
becomes 

1=0 u=0 

in which the root of unity is x'^C'^'' = C|f C''"' = = (2^"^"'^^'^ ■ Now, we use the 

fact that p2 = —1 mod 6, so that re{p2 + 1) = tp2 where P2 = {p2 + l)/<^- The above sum is now 

00 P2 — 1 00 

Y,Ck,iW'^(c'''''''-^'y =P2 Yl Ck,iw' = P2W-^'P'^ Yl (^k,jp2-ktp',q' ■ 

i=0 u=Q i=—ktp'2 mod p2 jl^^ip^/P'^ 

leading to the result. □ 

Proposition 4.11 In case pi = P2 = P, for all k ^ 0, 

sc,k e ( q-'/''^X'm] = Q'^'^'m], 



r]{z) 

where all series are explicited in the proof. 



Proof: One uses = Cp in 



sc,k = f'/' r-^) E^(-r'^'' 

V viz) J frt v{z + i^ 



viz + i^/pY'' 
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= p 

where 



efe p-1 kK{u) i-eke{v) 

, ^ ''24 

' (1- 

^2 \ efe P^l 



efc/2 / Nefc -efc/24A-efci//24^_^r_^_|,2^ 



:fc/2 /^^-l/24 ^(y^)^ ^ y^g^^-^efc-^fcK(;y)^efc6>(iy)^.-efci^/24^^^^^iy^efc 



ry(z) 



Writing C{qy^ = X^i^o Cfc,i'7* (same remark on the sign of A;), the inner sum of the preceding relation 



IS now 



p— 1 oo oo p— 1 

^,(^)efc^..H^efc.H^_eW24^^^_^(^^.). = ^ ^ ,(^)e/c^/c.(.)^e..H ^^-e./24^.).^ (3) 

!/=l 1=0 i=0 v=\ 



Let's treat the case p = 2 first, with e = 8. We get 

For p odd, the root of unity in the inner sum of ([3]) is 

the exponent of C24p being 

p{-{p - if^i + pv (1 - f?) + (-3p + 2 + t>) /u - 3 + 3p) - z/. 
When p = 3 and e = 3, we find 

M-mn+?,u-?,vp?+l?,) _ ,^-16^l+4ly-4u^l'^+9^k _ ^i^\k 

S.24 — \^12 I — Im ■ S3 J ) 

leading to 

When A: is even, this boils down to 



.3)3fe/2 / -1/24^'^'' 



i=—k mod 3 i^—k mod 3 

X 3fc / oo \ 

^ V •i=-fcmod3 i=0 / 



When k is odd 
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and 

'-i+k , ^2(i+k) _ ] \ii + k = Q mod 3, 

from which 

\2 \ 3fc oo 



2 

^ eli^J U3 j - ^3 +^3 - I (_i)(*+fc)mod3^ otherwise, 



i+fc^O mod 3 



When p > 3, we get 



ek 



/^-(p-l)V+P'^{l-/^^)+(-3p+2+j;)/i-3+3pNefc/A-efc+24iNi. _ ( {v-y IJ.'^ - IJ.+'i)p+tiv~Z \ , ^-ek+2Ai\u 
\'^2i ) \^24p ) — [^24 j \^24p ) 

using = 1 mod 24. We simphfy this as 

Mu+3)-3Y^ /^-efc+24^^l/ _ /-efe(p-l) / .efe(p2_l)+24j 

Write p2 _ 1 = 2V to obtain 
When A; is even, this gives 



\ j+efcp'=0 mod p j+efcp'^0 mod p 



^ /V / / Y i+efcp'=0 mod p i=0 

When k is odd, remarking that e is always odd from Table O the sum is now 

But ^IZ\ e(i^)(Cp^'^^V = when i + eV = mod p since there are the same number of quadratic 
residues and non-quadratic residues modulo p. When i + ekp' ^ mod p, Cjp^^^^ is a primitive p-th 
root of unity. Remember that |8| Ch. 6] 



E E y> 

X residue x non residue V ^ 

Let 5 be a generator of (TLjpTL^* . If u is an integer, then 



a: residue x non residue • 
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When i + ekp' ^ mod p, we set 0(i + ekp') = u such that g"^ = i + ekp' mod p. Then 



^C,k 




When p = 1 mod 4, the first terms simphfy to (^^^^^ ^•'/^ = -y^rhen p = 3 mod 4, we get 

_^2+efc(p-l) ^ _^l+efc(p-l)/2 ^ ^_-^-^{3+ek(p-l)/2)/2 

As a last point, the dominant term of Sc,k is When p = 2 and e = 8, this is 2k/3, 

whereas re = 1/3; when p = 3, e = 3, we get k/2, whereas re = 1/2. For p > 3, we have e = 1 and 
we compare {p — 1)/12 and re = {p — l)/2 • {p — 1)/12. Looking at the valuation of 2 and 3, we 
deduce that re = t/6 and {p — 1)/12 = p' /S. □ 

4.7 Tables of equations for double ?7-quotients 

^[ml2\{F, G2) = -G2F^ + 208 + 31 G2 F'^ + + 16, 
$[tt)|3] (F, G3) = _ ^11 _ i?io ^ 27 G3 - 10557 - 162 G3 F'^ - 14076 F^ - 18 G3 F^ 

-9801 F^ + 163 Ga F^ + (486 - Gi)F2 _ 9 G3 F - 27. 
$[tt)3,7](F, G3) = F^^-Gs f3^-514F3°+21 G3 f29-12585F28-147G3 f2'^-25158 f26+322 G3 F^^ 
-5103 F^^ + 378 G3 F^^ + 80556 F^^ - 1638 G3 F^^ - 21994 F^" - 28136 F^^ + 1620 G3 F^'^ + 25650 F^^ 
-252 G3 F^^ - 3944 F^^ - 322 G3 F^^ - 14938 F^^ _^ 22 G3 F^^ - (G^ - 2940)F1° - 10 G3 F^ + 1953 F^ 
+G3 F"^ - 462 F^ + 7 G3 F^ + 15 F^ - G3 F^ - 10 F^ + 1. 

5 Conclusion 

We have studied modular equations involving 72 and 73 for double r/-quotients. As a result, more 
compact modular equations can be stored and used, with application to the SEA algorithm (see 
for instance |T0]). or CM computations, as motivated for instance by [13] (see |llj). 

It seems natural to conjecture that more general functions can exhibit the same properties. 
Experiments can be conducted on Newman functions, using for instance the resultant approach, 
leading to new instances of the theorems. This will described in another article. 
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